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Ultrasound velocity measurements of the orbitally frustrated spinel GeCo2 04 reveal colorful and 
unusual elastic anomalies within the antiferromagnetic phase. Temperature dependence of shear 
moduli exhibits a minimum within the antiferromagnetic phase, suggesting the coupling of shear 
acoustic phonons to molecular spin-orbit excitations. Magnetic-field dependence of elastic moduli 
exhibits diplike anomalies, being interpreted as magnetic-field-induced metamagnetic and structural 
transitions. These elastic anomalies suggest that the survival of geometrical frustration, and the 
interplay of spin, orbital, and lattice degrees of freedom evoke a rich set of phenomena in the 
antiferromagnetic phase. 

PACS numbers: 71.70.Ch, 71.70.Gm, 75.25.Dk, 75.50.Ee 



The concept of geometrical frustration provides an in- 
triguing playground for condensed matter physics. So 
far, the majority of studies have been devoted to the spin 
degrees of freedom^. Here, we would like to focus on the 
orbital sector as the leading role. An orbital system is in- 
herently frustrated even on a simple square lattice: when 
orbitals (directions of the electron cloud) are arranged to 
gain bond energy for one direction, this configuration is 
not fully favorable for other bonds^Tii. Furthermore, the 
orbital frustration is expected to be enhanced on a ge- 
ometrically frustrated lattice, where the orbitals as well 
as spins cannot form antiferro-type orderi^*^. 

Cubic spinels AB2O4 with magnetic B ions have at- 
tracted much interest in light of the geometrical frus- 
tration which is inherent in the B-site sublattice of cor- 
ner sharing tetrahedra. The spinel cobaltite GeCo204 
consists of magnetic Co^+ {3d^) on the octahedral B 
sites and non-magnetic Ge^+ on the tetrahedral A sites. 
The positive Weiss temperature, Qw = +81.0 K, indi- 
cates the dominant contribution of ferromagnetic (FM) 
interactions^. However, an antiferromagnetic (AF) long- 
range order appears below Tn = 21.0 K, in coincidence 
with a cubic-to-tetragonal structural elongation with c/a 
= 1.001—"—. The AF order is composed of a dominant 
component described by a trigonal propagation vector 
Qii = (1/2, 1/2, 1/2) and an additional one described by 
a tetragonal Qj = (1, 0, 0)/(0, 1, 0), where the magnitude 
of Qi is much weaker than that of Qjj^^^^. 

GeCo204 is a promising candidate for the frustrated 
system with orbital degrees of freedom^^'^^. Ultra- 
sound velocity measurements suggested that the struc- 
tural transition at Tn in GeCo204 is ascribed to the re- 
lease of frustration rather than the Jahn- Teller eSeci^. 
An octahedral-site Co^+ is a well-known spin-orbit (SO)- 
coupling active ion with an unquenched orbital angu- 
lar momentum L. The electronic ground state is nor- 



mally described by effective total angular momentum 
Jeff = 1/2 doublet with L = 1 based on the triply- 
degenerate t2g orbitals and S = 3/ 2-^^1-^^ . In GeCo204, 
the Jeff = 1/2 ground doublet was confirmed by spe- 
cific heat measurements and inelastic neutron scattering 
experimentsi^ii^. 

An inelastic neutron scattering study discovered that 
a strong frustration effect (molecular spin excitations) 
still survives even below Tn in the typical spin-frustrated 
spinel chromite MgCr204— . Very recently, it was also re- 
ported that GeCo204 exhibits molecular SO excitations, 
not only above but also below Tat^. Therefore, novel 
phenomena relevant to the geometrical frustration are 
expected to be hidden in the AF state in GeCo204. 

In this paper, we present ultrasound velocity measure- 
ments in GeCo204. The modified sound dispersions by 
the magneto-elastic coupling allow one to extract detailed 
information about the interplay of lattice, spin, and or- 
bital degrees of freedom. We find colorful and unusual 
elastic anomalies in GeCo204 in the AF state. 

Ultrasound velocities were measured in single crystals 
of GeCo204 prepared by the floating zone methodic, 
where the phase comparison technique was utilized with 
longitudinal and transverse sound waves at a frequency 
of 30 MHz. The ultrasound waves were generated and 
detected by LiNbOa transducers glued on the parallel 
mirror surfaces of the crystal. We measured sound veloc- 
ities in all the symmetrically independent elastic modes 
in the cubic crystal: compression modulus Cu {Aig sym- 
metry), tetragonal shear modulus = Ct {Eg 
symmetry), and trigonal shear modulus C44 {T2g sym- 
metry). 

Figures 1(a) and (b) depict the transverse sound ve- 
locity Vt in C44 and Ct with magnetic field i^||[110] as 
functions of temperature (T), respectively. The theoret- 
ical extension of the experimental data in the paramag- 
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FIG. 1: (Color online). T dependence of vt and vl with 
if II [110]. (a) C44, (b) ^^^^-^^^^ = Ct, and (c) Cn. The 
dotted curve indicates the background Cij in each modulu^^^. 
The absolute value of the sound velocity at 2 K is t'?^ = 4030 
m/s for C744, = 3300 m/s for Ct, and vl = 6940 m/s for 
Cii, respectively. 



netic (PM) region of 50 K < T <150 K down to T ^ 0, 
the background C^j , is also shown as the dotted curvei^. 
As reported in our previous study^^, C44(T) and CtiT) 
with H = exhibit, in addition to a diplike anomaly at 
Tat = 20.6 K, nonmonotonic T dependence in T < Tat 
with a characteristic minimum at ~6 K^^. Furthermore, 
as a newly discovered feature, Ct{T) at 2 T exhibits the 
steep softening below ~6 K with decreasing T. 

We explain that the minimum in C44(T) and Ct{T) 
with if = 0, and the softening in Ct{T) at 2 T arise 
from magnetic origins, not from domain- wall effect. In 
the magnetically ordered state, magnetostriction can give 
rise to the domain- wall stress which is detected as loss of 
elasticit}4i^. This stress is generally removed by applying 
if, which is observed as recovery of elasticity with H. For 
GeCo204, the domain- wall stress effect on the volume- 
strain mode Cii should be stronger than the volume- 
conserving modes C44 and Ct. Fig. 1(c) depicts T de- 
pendence of the longitudinal sound velocity vl in Cu. 
Below Tat, Cii(T) with H = loses the elasticity com- 
pared to 3 T. Thus Cii(T) with H = should be affected 
by the domain- wall stress. However, Cii(T) with H = 
is almost constant below Tat, indicating the T indepen- 
dence of the domain- wall stress effect. 

Figures 2(a) and (b) depict the relative change of vt in 
C44 and Ct as functions of ii'||[110], respectively. Cm{H) 
and Ct{H) exhibit several distinct anomalies within the 
AF phase at 2.0 K and 17.5 K, in contrast to the absence 
in the PM phase at 30 K. At 17.5 K, a diplike anomaly is 
seen in Cm{H) and Ct{H) at ijlqHn = 6.6 T as marked 
by the open arrows, which coincides with the AF-to-PM 
transition observed in the magnetization as well as the 
specific heat measurement92.5£^. Below Hn at 17.5 K, 
C44{H) exhibits steep hardening with increasing ii" in '-«^1 
T< /JjqH <~3 T, in contrast to the weak H dependence 
in CtiH). At 2.0 K, if at shifts to over 7 T, and further 
anomalies evolve as marked by the filled/open circles. 

Although the isothermal magnetization measurements 
at 1.5 K with if 1 1 [110] reported a step-like anomaly at ~4 
T suggesting the if -induced transition^, no Cij{H) at 2 
K in Figs. 2(a)-(c) exhibits anomaly at ~4 T. Very re- 



cently, neutron scattering experiments under if identified 
this ~4 T anomaly as a spin reorientation transition in 
Qjj^. The absence of anomaly in Cij{H) suggests that, 
in contrast to the magnetic susceptibility —d^F/dH'^, 
the strain susceptibility d'^F/de^ is insensitive to the if- 
induced transition at ~4 T in Qn. In fact, the trigonal 
Qii is inconsistent with the tetragonal lattice distortion 
in symmetry, suggesting weakness of the coupling be- 
tween Qii and lattice. 

Assuming that the domain- wall stress effect in di(if) 
at 2.0 K is expressed by the linear fit in < iiqH < 0.5 T 
as indicated by the dashed line in Fig. 2(c), Figs. 2(d)- 
(f) depict C44(i^), Ct{H), and Cii(ff) at 2.0 K with the 
subtraction of the domain- wall stress effect, respectively. 
C44 (if ) exhibits the hardening with increasing if in ~0.5 
T< fioH <-3 T. Referring to Cm{T) in Fig. 1(a), this 
hardening in Cm{H) corresponds to the recovery towards 
the background elasticity with if, C44 — C44 = AC44 
0. Thus the most probable origin for the hardening in 
C44 (if) is the suppression of the magneto-elastic cou- 
pling with increasing if. In contrast, Ct(if) exhibits the 
softening with increasing if in < /j^qH <~6 T as indi- 
cated by the dotted line in Fig. 2(e), which corresponds 
to the increase of — Ct = ACt with increasing if. 
Thus the softening in Ct{H) indicates that, in contrast 
to C44{H)^ the magneto-elastic coupling is enhanced with 
increasing if. Such elastic-mode-dependent if variations 




FIG. 2: (Color online), (a)-(c) The relative shift of vt and 
VL as functions of if ||[110]. (a) C44, (b) Ct, and (c) Cu. The 
curves are shifted for clarity. The open arrows in (a) and (b) 
indicate the AF-to-PM transition at 17.5 K. The dashed line 
in (c) is a linear fit to Cii(if) at 2.0 K in < /xoif < 0.5 T. 
(d)-(f) The data at 2.0 K in (a)-(c) after the subtraction of 
the domain- wall stress effect. The dotted line in (e) is a guide 
to the eye indicating the softening in Ct(if). Inset shows the 
magnetization curve at 2.0 K. The filled/open circles in (a)- 
(f) indicate the elastic anomalies with/without the magnetic 
hysteresis. 



3 



of the elastic constants, the hardening in C4^4^{H) and 
the softening in Ct{H)^ are observed most Hkely as the 
results of the coupling of the sound waves to the spin- 
wave excitations^^. At = 0, referring to Figs. 1(a) 
and (b), AC44/C44 is much larger than ACt/Ct, indicat- 
ing that the spin- wave excitations affect C44 much more 
strongly than Cf. The downturn softening in C44{H) and 
Cii{H) in ~3 T< jnoH <7 T would be the precursor to 
the AF-to-PM transition at the higher H. 

In Figs. 2(d)-(f), another salient feature is the diplike 
anomalies marked by the filled/open circles, indicative of 
the i^-induced transitions. These anomalies evolve in the 
lower T in the AF state, and are elastic-mode selective 
with the occurrence at ^1 T in all the elastic moduli, at 
~2 T in Ct and Cn, and at ~3 T in Cu. Here, we regard 
the shoulder structure in C44{H) at ~1 T as a diplike 
anomaly because this structure should be formed by the 
superposition of the hardening in ~0.5 T< /j^qH <~3 T 
on the diplike feature. The diplike anomaly in Cij{H) 
within the magnetically ordered state is normally ob- 
served at the metamagnetic transition^^. The inset to 
Figs. 2(d)-(f) shows the magnetization curve at 2.0 K 
with i7||[110], where the change of slope and the hys- 
teretic feature are seen in < jj^qH <~1 T. Thus the 
elastic anomalies at ~1 T in Figs. 2(d)-(f) marked by 
the filled circles should be attributed to the metamag- 
netic transition. In contrast, the elastic anomalies at ~2 
T and ~3 T in Figs. 2(e) and (f) marked by the open cir- 
cles are attributed to other i^-induced transitions with- 
out anomaly in the magnetization curve. 

From now on, let us discuss in detail the origins of the 
two types of unusual elastic anomalies observed within 
the AF phase: a characteristic minimum in C44(T) and 
Ct{T) with = 0, and diplike anomalies in C44(i^), 
Ct{H), and Cii{H). First, we discuss the origin of a min- 
imum at ~6 K in C44(T) and Ct(T) with H = shown 
in Figs. 1(a) and (b). According to the specific heat 
measurements in GeCo204, there is no additional phase 
transition within the AF phase with H = 0^. Thus 
the minimum in C44(T) and CtiT) with H = should 
originate from the coupling of the sound waves (acoustic 
phonons) to the magnetic excitations, not to the static 
order. As the possible magnetic excitations, we can first 
consider the spin- wave excitations. In GeCo204, as de- 
scribed above in conjunction with Figs. 2(d) and (e), the 
coupling of the sound wave to the spin-wave excitations 
should be suppressed in C44 but enhanced in Ct with 
increasing H. However, as exemplified by the data at 
3 T in Figs. 1(a) and (b), the minimum disappears by 
the application of > 2 T in not only Cm{T) but also 
Ct{T)^. Thus we rule out the spin- wave excitations as 
the origin of the minimum in C44(T) and CtiT) with H 
= 0. We now note that the inelastic neutron scatter- 
ing study in GeCo204 discovered the appearance of an 
energy gap A ~ 3 meV as the low-energy magnetic ex- 
citations in the AF statei^ii^. And, interestingly, this 
3 meV-gapped mode was characterized as the molecu- 
lar SO excitations which coexist with the AF order—. 



Therefore, the minimum in C44{T) and Ct{T) with H 
= at T < Tn should stem from the molecular SO ex- 
citations. Noteworthily, this conclusion means that the 
elastic anomalies due to the molecular magnetic excita- 
tions (frustration effect) have been observed, as far as we 
know, for the first time. 

Here, we would like to give a quantitative discussion of 
the molecular SO excitations A using CtiT) with H = 
in Fig. 1(b), shown again in Fig. 3(a) (symbols). The 
inelastic neutron scattering study in GeCo204 suggested 
that, in the ground state, the nonmagnetic singlet pairs 
of the FM Co^+ tetrahedrons, " di-tetramers" shown in 
the inset to Fig. 3(a), coexist with the AF order Here, 
A ~ 3 meV (~35 K) corresponds to the singlet^triplet 
excitation energy of the di-tetramers^^. A characteristic 
minimum in CtiT) with H = due to such short-range 
magnetic excitations is analogous to the one in the low- 
dimensional spin dimer systems like SrCu2(B03)2^^. The 
contribution of the di-tetramers to the elastic constant is 
written as, 
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Xs 



(1) 



with the density of di-tetramers, C = \dA/de\ the 
coupling constant to a single di-tetramer which measures 
the strain dependence of the singlet- triplet gap, K the 
inter-di-tetramer interaction, and Xs = G~'^{d'^F / de^) 
the strain susceptibility of a single di-tetramer. Since 
A is expected to evolve in the lower T from the un- 
gapped state at T^^^^^ we assume the T-dependent 
gap, A(T) = A(l - {T/Tn)""). Here, the value of A = 
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FIG. 3: (Color online), (a) T dependence of Ct with H — 
from Fig. 1(b) (symbols). The solid curve is a fit to Eq. (1) in 
assumption of the molecular SO excitations A = 35 K. Inset 
picture shows a di-tetramer in the B-siie Co^^ sublattice. 
The fit parameters are summarized in the right table, (b) 
if-induced metamagnetic and structural transitions in the Qi 
order schematically depicted in a Co^^ tetrahedron. The filled 
and open circles depict up and down spins. 
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35 K is fixed as the first approximation. In Fig. 3(a), it 
is evident that a fit to Eq. (1) shown as the sohd curve 
excehently reproduces a minimum at ~6 K in CtiT) sup- 
porting the presence of the molecular SO excitations. 
The large fitted value oi G = 1334 K indicates that the 
dimerized FM tetrahedrons (di-tetramers) in the tetrag- 
onal AF phase strongly couple to the lattice deformation. 
On the other hand, the absence of anomaly in CtiT) in 
the PM phaseii indicates G for the gapless molec- 
ular state of the undimerized FM tetrahedrons in the 
cubic PM phasei^. These features imply the correlation 
between the lattice structure and the SO molecular state. 

Second, we discuss the origins of the diplike anomalies 
in Caa{H), Ct{H), and Cii{H) at 2 K shown in Figs. 
2(d)-(f). For a sound wave with polarization u and prop- 
agation /c, the magneto-elastic coupling acting on the ex- 
change interactions is written a^, 

He.s = ■ ")(^ • '^)(^^ • ^.Oe^C^ ^--*). (2) 

i 

Here S = Ri — Rj is the distance between two magnetic 
ions. On the basis of this formula, the metamagnetic 
transition at -1 T observed in C44 (^||[001], ii||[110]), Ct 
(A:||[110], ix||[110]), and Cu {k\\u\\[100]) should be dom- 
inated by the exchange interactions in <110> direction 
because the exchange interactions in <100> and <111> 
directions become inactive in C44 and Ct, respectively. 
Thus the metamagnetic transition at ~1 T should oc- 
cur in the Qj order which is composed of the antiferro- 
magnetic spin chains along [110] and [110] directions as 
shown in Fig. 3(b)(i)^^. Note that this metamagnetic 
transition leads to the disappearance of the Qj order, as 
shown in Fig. 3(b) (ii). Therefore, it is expected that the 
tetragonal-to-cubic structural transition coincides with 
the metamagnetic transition at ~1 T, as shown in Figs. 
3(b) (i) and (ii). 

The iJ-induced transition at ~2 T is characterized by 
the anomaly in Ct{H) as shown in Fig. 2(e) (open cir- 
cle), the evolution of which corresponds to the softening 



in Ct{T) at 2 T below -6 K, as shown in Fig. 1(b). The 
absence of anomaly at ~2 T in the isothermal magneti- 
zation, as shown in the inset to Figs. 2(d)-(f), rules out 
the possible metamagnetic transition. Thus the anomaly 
in the tetragonal Ct at ~2 T strongly suggests the oc- 
currence of the cubic-to-tetragonal structural transition, 
as shown in Figs. 3(b) (ii) and (iii). The anomaly at ~3 
T only in Cii{H) {Aig symmetry) suggests, with the ab- 
sence of anomaly in the isothermal magnetization^, the 
occurrence of the isostructural non-metamagnetic tran- 
sition. The nature of this transition is uncovered at this 
stage. Taking into account the SO-coupling active Co^+, 
it is expected that the i^-induced structural transition at 
~2 T coincides with the orbital state transition driven by 
H and SO coupling. Detailed studies by applying the mi- 
croscopic probes such as NMR, neutron scattering, and 
resonant X-ray scattering would uncover the exotic or- 
bital phenomena. 

As mentioned above in conjunction with Figs. 1(a) 
and (b), the minimum in Ct(T) and C44(T) with H = 
disappears in the data with i7 > 2 T, indicating the 
decoupling of the shear acoustic phonons from the molec- 
ular SO excitations. This i^-induced decoupling might 
be driven by the JY-induced structural transition at ^2 
T. Here, the transformation or the disappearance of the 
SO molecules is expected to coincide with the i7-induced 
structural transition. 

In summary, ultrasound velocity measurements of 
GeCo204 reveal unusual elastic anomalies within the AF 
phase due to the molecular SO excitations and the H- 
induced metamagnetic and structural transitions. The 
present results suggest that the survival of geometrical 
frustration, and the interplay of spin, orbital, and lat- 
tice degrees of freedom evoke a rich set of phenomena in 
the AF phase which deserves further experimental and 
theoretical studies. 
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